UNIVERSITY OF ILLINOIS 
DIGITAL COMPUTER 



LIBRARY ROUTINE H6 - 86 



TITLE 

TYPE 

NUMBER OF WORDS 

TEMPORARY STORAGE 



ACCURACY 



DURATION 



READ AROUND 



.PRESET PARAMETERS 



S3 


OOF 00 aF 


S^ 


OOF 00 bF 


35 


OOF 00 cF 


s6 


OOF 00 d^J 


ST 


OOF 00 &T 


S8 


OOF 00 sF 


S9 


OOF QO(a+n)F 


SK 


OOF 00 (2b )F 


ss 


OOF 00(c+n)F 



Minimization of a Function of n Variables "by Treating One 
Variable at a Time (DOI or SADOl) 
Closed 
k6 

n+1 "words at S3 

1 word at Sk 

n+1 words at S5 

Depends upon the condition of the function. In general 
the ultimate error in the position of the minimum will 
be less than the last mesh size used, (See description) 
Depends upon the condition of the function. In general 
the duration will only be a trivial amount longer than 
the time to compute the function f (x^x^, .. «* x n ) the 
great number of times necessary using this crude (see 
description) approach. 

Depends upon the condition of the function and the routine 
used to compute it. The read around will be a little less 
than that of the function subroutine when used continuously. 
S3 - SS 

Where a 5 a+l f ..*, a+n are the addresses of f(x ,*.*,x n ) 
x-, Xgjr-.-jpc on entry and exitf 
b is the address of b, the current mesh sizef 
c, c+1, ... y c+n are working spaces^ 
^C is a factor by which the mesh size is decreased^ 
e is an end constant such that only mesh sizes larger 
than e will be used§ 

s is the address of a closed subroutine which takes x^ 
from c+l(lS5) jr x 2 from c+2 (2S5 )***», x n from c+n(nS5) 
and places f (x-y>L, ...^) in R^ 
n is the number of variables . 



ENTRY 



x , x p> •>*, x the initial values of x 1 j,x p - t * . « ,x 
are to "be loaded in 1S3, 2S3, •••,nS3- S , the initial 
value of 5 * is to be in R. ,■ and the main program should 
contains 



p 


50 


rF 




50 


pF 


p+1 


26 


:.QF 



■where q is the address of this routine and r is the address 
to the left hand side of which control will be transferred 
before every decrease in mesh size. At this address the 
programmer may place a routine to assess or print inter- 
mediate results * This routine should return control to the 
left hand side of (q+32)» This can be done automatically 
since when control is transferred to r(L.Ho) the appropriate 
return address will be in the right hand address position 
of R p and can be utilized by S5, ^2 into a 26 at the end 
of the programmer's interlude routine. The best current 
values of f(x , ...pc ), x ,.. ».,x are always in S3, 1S3*»*«.» 
nS3- If it is not desired to leave this code, set 
r==Ji+32* The routine is finally left with the best values 
of f(x 1 ^x 2 , «. i x n ), x x , Xg,...,*^ in S3* 1S3, 2S3,-*»,nS3« 
METHOD 

1- The routine examines f(x_. #•*.•>:£ ) at x, .+ 5 going in the 

decreasing direction until f is minimized as a function 
of x r 

2- Step 1 is repeated for x_ ,*,.** *«>x_' in sequence. 

3- We then return to x $ then x p , etc- until no improvement 

would be obtained for x., + 5* x« + 5 » . »** x .+ 6 » 

1 — * 2 — ' * J n — 

%- Control is then transferred to the interlude. 

5- Upon returning to the routine, the current mesh size, 5., 

is replaced by 6. =1^8. where o( is the number entered 
in S6 during read-in o Steps 1-4 are then repeated . 
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6- Steps 1-6 will be repeated (N+l) times where N is the largest 

integer for which «^ 8 >e. € is to he entered in SJ during 
read-in. 

NOTE A function,,,, poorly conditioned (e.g.) in having a very small 

gradient with respect to some argument, can deceive a method 
of steepest descent routine in that this coordinate of the 
minimum will be very poorly found « If the programmer suspects 
that such a condition exists^ Code 85 can be used for roughing 
and this code, a brute force approach to the minimum by varying 
one argument at a time may be used for finishing- All 
parameters (except SK and SN), entry* and contents of storage 
location of this code and Code 85 are identical, so that they 
may be used at the same time in the machine for the above 
purpose * 



RtS 7/22/59 



BAT E 10/2/56 RT.g 9/12/58 
CODED B Y Jo N* Snyder 
APPROVED BY J* P* Mash 
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LOCATION 


ORDER 


HOTES MGE 1 





QpK (H6) 

ko sff 




Store 5 




S5 32L 




Link in 


1 


46 15L 
L4 3% 




Plant interlude address 


2 


42 3^L 




Plant link address 




L5 38L 




Set transfer 


3 


4o 4l 

50 211 




Orders for initial entry 


4 


L5 F 


From 43 


x, , x_ j * * « yX. to 1S5, * » * ,nS5 




40 F 


By 3,6 




5 


L5 4l 




Advance transfer orders 




L4 35L 




- 


6 


40 4l 








L5 3% 






T 


LO 4l 




■Test for i = n 




26 43L 




_ 


8 


L5 40L 


From 31,34 " 






42 13L 


and 45 


Set routine to process x ± 


9 


42 24L 
46 15L 






10 


46 25L 








40 19L 


- 


1 


11 


L5 13L 




" Close switch for complete failure -with 




46 31L 


- 


given mesh size 


12 


L5 3L 


From 30 


Tciose switch for failure in one direction 




42 171 






13 


50 F 


By 1 


Interlude address 


■ 


L5 F 


By 8,27 
f rem 24 


Form x. + 5 

1 — 


14 

1 


LO S4 
50 F 


By 22 




15 


40 F ' 


By 9,28 


f (,..«,x. + 5 ,.,*«,) to S5 




50 15L 




- 


.16 


26 S8 








40 S5 






IT 


LO S3 




f - f (hefore) 




32 F 


By 12,23 


Directional failure switch 



86 



LOCATION 


ORDER 


NOTES VkBR P 


18 


L5 S5 






Replace Old values by new ones 




40 S3 


. 




r ' 


19 


L5 F' 


By 10, 26 






. 


40F 








20 


L5 34L 






Open mesh size 

""" '■"■ 




46 31L 






failure switch 


21 


26 23L 








L5 42L 


From I? 




•■ 1 
Replace + "by - or - "by + 


22 


LO 14L 






. 




40 14L 






■ ' 1 


23 


L5 3H 


From 21 




Open directional failure switch " 




42 17L 


1 




. 


24 


22 13L . 


From 17 


■ 




L5F 


By 9, 27 


Replace falling values 


25 


40 F 


By 10, 28 


■' 




L5 19L 








26 


L4 35L 






» ■ 




40 19L 






Advance addresses for treatment of 


27 


42 15L 








28 


42 24L 
46 15L 
46 25L 




> 


X i+1 


29 


L5 4lL 
LO 19L 






Test for i * n 


50 


56 12L 










50 L 




Set up link address 


31 


26 F 


By 11, 20 


Mesh size failure switch 




00 24L 




■ 


52 


50 S4 
13 36L 


■ 

!■ 1 




Form S . , =»^ S . and test against € 


33 


40 S4 
LO 57L 






' 


3^ 


36 8L 










22 F 


By 2 


Link Address 


35 


00 IF 




■ 




00 IF 




Advancer 



86 



LOCATION 


ORDER 


NOTES PAG$ 3 _ / 


36 


00 F. 
00 S6 




■ O^ 


37 


00 F 








00 S7 




* € 


38 


L5 1S3 
40 1S5 




Starting constant 


39 


L5 S9 
40 SS 




End constant 


4o 


L5 1S5 
40 1S3 




Starting constant 


4i 


L5 SS 








40 S9 


• 


End constant 


42 


F4 SK 








KD F 




(+) + (-) 


43 


36 4l 
50 43L 






f(x°,...,x£) to S3 


44 


26 S8 
40 S3 






— 


45 


26 8l 










00 F 


■ 


■ 



